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SUMMARY 


Measurements have been made of the mean total-head and temperature 
fields in a round turbulent jet with various initial temperatures. The 
results show that the jet spreads more rapidly as its density becomes 
lower than that of the receiving medium, even when the difference is not 
sufficiently great to cause measurable deviations from the constant- 
density, dimensionless, dynamic— pressure profile function. Rough 
analytical considerations have given the same relative spread. 

The effective "turbulent Prandtl number" for a section of the fully 
developed Jet was found to be equal to the true (laminar) Prandtl number 
within the accuracy of measurement. 

Measurements of turbulence level (u', v'), temperature fluctuation 
level d * , and temperature— velocity correlation du permit a comparison 
of their relative magnitudes. 

Direct measurements have been made of the double correlations uv 
and dv across a section of the fully developed jet, and the shear- 
stress and heat— transfer distributions have been computed therefrom. 
Finally, these last-mentioned measurements have permitted a determina- 
tion of the distribution of turbulent Prandtl number across the Jet, 
and these values agree quite well on the average with the effective 
value computed from mean velocity and temperature alone. 


INTRODUCTION 


The present work is concerned with two particular problems on the 
flow of round, low— velocity, turbulent Jets: (a) The effect of mean- 

density differences upon the rate of spread of the jet, which can be 
examined through the average velocity and density fields; (b) the rela- 
tive rates of transfer of heat and momentum in a heated jet , which must 
involve a study of fluctuations in velocity and temperature as well as 
average fields. 

Problem (a) has been investigated experimentally by Pabst (refer- 
ences 1 and 2), Von Bohl (reference 2), and others. Some British results 
consist of a few rough measurements made in the course of a different 
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investigation and, unfortunately, the German, reports are available only 
in practically illegible form, so that a thorough reading has not been 
feasible. These references arrived when the present investigation was • 
in progress and substantiated the result, discussed in detail in the 
body of the paper, that reduction in density of jet relative to receiving 
medium increases the rate of spread. 

Theoretical ana lysis of the two-dimensional, variable-density jet 
has been attempted by Hu (reference 3)* using the momentum-transfer 
theory. However, he has begun by assuming the constant— density charac- 
teristics of linear Jet spread and inverse parabolic decrease in axial 
velocity and therefore has not actually solved the variable— density 
problem. 

Goff and Coogan (reference 4 ) and Abramovich (reference 5) have 
solved the two-dimensional, variable— density, single-mixing zone, by 
momentum— transfer and vorticity— transfer theories, respectively. The 
former analysis involves the drastic assumption of a discontinuous 
density field. 

In general, however, the shortcomings of the theory of mixing length 
and the more recent theory of constant exchange coefficient, discussed 
briefly by Liepmann and Laufer (reference 6 ), render such analyses useful 
perhaps only for qualitative results. In fact, as illustrated in refer- 
ence 6 , equally useful results in shear flow problems are obtainable with 
appreciably less difficulty by the use of the integrated equations of 
motion, with reasonable guesses for the shape of velocity and/or shear 
profiles. 

A related method has been used by Ribner (reference 7 ) , who has given 
an approximate solution for the variable-density round jet in a moving 
medium by a rather empirical generalization of the method used by Squire 
and Trouncer (reference 8 ) in constant-density Jets. However, this method 
employs the shear assumption of the momentumr-transfer theory, so that the 
result, though a saving in labor, cannot be expected to give greater 
accuracy than the complete momentum— transfer analysis. 

In the present report, a very brief and approximate analysis is 
given for the relative jet spread as a function of density ratio, 
obtainable without a complete solution of the integral problem. 

The problem just noted is different from that of a jet of high 
subsonic velocity, in which the temperature differences arise from 
frictional heating. However, since the maximum rate of dissipation'of 
turbulent energy (into heat) in a jet takes place on the jet axis (refer- 
ence 9 ), the overall behavior of a fully developed, turbulent, high- 
speed, subsonic jet may not be very much different from a lower speed 
heated Jet, if the high-speed jet starts out with the same temperature 
as the receiving medium. There appear to be no data available on 
turbulent jets with Mach number approaching unity. Abramovich (refer- 
ence 5) has also applied the vorticity— transfer theory to an approximate 
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solution of the high-speed (subsonic), plane, s ingle— mixing region with 
the same stagnation temperature for moving and stationary mediums. 

The present investigation on problem (a) has been undertaken 
primarily to determine experimentally the effect of density difference 
upon the jet spread. Another matter of interest is the possible deviation 
from simple geometrical similarity for large density differences. 

Somewhat more work has been previously done on the simple phases of 
problem (b), the relative rates of heat and momentum transfer in a round 
turbulent jet. The first systematic measurements of velocity and temper- 
ature distribution were apparently made by Ruden (reference 10 ), who 
found that the temperature distribution is appreciably broader than the 
velocity distribution in the fully developed region. A more detailed 
investigation of the kinematic quantities, including turbulence level 
near the axis, was carried out by Euethe (reference 11 ). The results of 
these two investigations have been verified and extended by one of the 
present writers (reference 12). 

Since there exists no satisfactory theory of turbulent shear flow, 
naturally there is none for heat transfer in turbulent shear flow. Most 
of the theoretical work has been done on channel and boundary— layer flows, 
with the heat transfer set up in direct analogy to the momentum transfer. 

The empirical factor connecting the two is, of course, exactly the quantity 
that a satisfactory theory must predict. 

For turbulent jets, momentum— transfer considerations lead to identical 
curves for velocity and temperature distributions (reference 13 ) , whereas 
the modified vorticity— transfer theory leads to the qual.it at ively correct 
result (reference l4 ) that heat diffuses more rapidly than momentum. 
Quantitatively, however, this result is considerably in error (reference 12). 

Clearly, no satisfactory investigation, either experimental or 
theoretical, of a turbulent flow phenomenon can be made without considera- 
tion of the fluctuations as well as the mean distributions. In the round 
turbulent jet, fairly detailed measurements have been made of the fluctua- 
tion kinematic quantities (reference 12), but no previous investigations 
appear to have included the significant quantities involving temperature 
fluctuations. Since a hot-wire technique for the simultaneous measure- 
ment of temperature and velocity fluctuations has recently been developed 
(references 15 and l6), it is now possible to study turbulent heat 
transfer from the points of view of both mean and fluctuating variables. 

At the present time, a weakness of the new measuring technique is 
the uncertainty of the detailed physical form of King*s equation (refer- 
ence 17) for the heat loss from a fine wire to a flowing fluid. A brief 
discussion is given in appendix A of reference 1 6 . The results presented 
in the present report indicate that, at least for variable temperature 
applications, King*s equation may be quantitatively satisfactory. 

The first measurements of temperature fluctuation level were made 
in 1946, in connection with the preliminary phases of the present 
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investigation. However, deviations from rectangular ity in the initial jet 
temperature distribution prompted the construction of a completely new hot- 
jet unit, used in all measurements presented herein, and the preliminary 
measurements have not been published. 

Measurements of mean total-head and temperature distributions at 
various jet temperatures have been made in order to find out if there is 
any appreciable change in the relative rates of heat and momentum 
transfer with absolute temperature, over a moderate range. 

Finally, it should be noted that investigations of heat transfer in 
turbulent shear flow are important not only for the immediate results but 
also as a means of studying the turbulent motion itself. 

This investigation was conducted at the California Institute of 
Technology under the sponsorship and with the financial assistance of the 
National Advisory Committee for Aeronautics. The authors would like to 
acknowledge the assistance of Dorothy Kerns, Beverly Cottingham, 

Sally Rubsamen, and Betsy Barnhart in computing the results and drawing 
the final figures for this report. 
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SYMBOLS 

diameter of orifice (l inch) 
axial distance from orifice 

radial distance from jet axis 

axial component of mean velocity 

radial component of mean velocity 

maximum U at a section on jet axis 

maximum U in the jet (i.e., in the potential cone) 

dynamic pressure ^ ptF’y 

maximum dynamic pressure at a section 

maximum dynamic pressure in the jet 

mean absolute temperature at a point in the flow 

absolute temperature of the receiving medium 
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absolute initial jet temperature 
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maximum 9 at a section 
maximum 9 In the jet (T 0 — T r ) 

axial component of instantaneous velocity fluctuation 
radial component of instantaneous velocity fluctuation 
instantaneous temperature fluctuation 
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P 
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P* 


total pressure 
air density 

mean density at a point in the flow (at T a ) 
density fluctuation 


Pmin 


Pc 

P 

V 

k 


minimum p at a section on the axis 

minimum p in the jet (at T Q ) 

density of receiving medium (at T r ) 

viscosity coefficient of air 

kinematic viscosity (p/p) 

thermal conductivity of air 

specific heat of air at constant pressure 
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Prandtl number 
shear stress 

heat transfer per unit area 
turbulent exchange coefficient 
turbulent heat— transfer coefficient 

"turbulent Prandtl number" 

average or effective cr^ across a section of the Jet 



A 


"momentum diameter" of Jet 



A e 


A at a section in a constant-density Jet 
"thermal diameter" of hot Jet 
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double correlation coefficients at a point 




UV 

U , Y f 
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instantaneous hotwire voltage fluctuation 

sensitivity of skew wire to u, and v fluctuations, 
respectively 
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EQUIPMENT 

Aerodynamic Equipment 


The 1— inch hot— jet unit is shown schematically in figure 1. The 
three— stage axial blower (reference 18) is driven by a 2-horsepower 
variable— frequency motor, operated at a fraction of rating. Heat is 
added through two double banks of coils of No. 1 6 Nichrome wire. As 
can be seen in the sketch, a good part of the heated air is ducted 
around the outside of the jet— air pip© in order to maintain a flat 
initial temperature distribution in the jet. A blower is used to help 
the air through this secondary heating annulus. Since this auxiliary 
air is not ducted back into the main blower intake, heat is discharged 
from the annulus at an appreciable rate and is exhausted outside of the 
building to minimize room-temperature rise. 

The relatively high velocity section between heaters and final 
pressure box promotes adequate mixing behind the grid to insure a 
uniform initial jet— temperature distribution. An earlier, simpler unit, 
mentioned in the INTRODUCTION, had the heating coils in the final (and 
only) pressure box; as a result, the Reynolds number of the coil wire 
and even that of an extra ceramic grid were too low to produce turbulent 
mixing. 

The complete unit is photographed in figure 2(a), and figure 2(b) 
is a close-up of the asbestos orifice plate and the traversing mechanism. 

Runs have been made at orifice velocities between 65 and 115 feet 
per second and at orifice— temperature differences between 0° C and 3^5° C 
The initial dynamic— pressure distribution (fig. 3) is effectively 
rectangular, and the initial temperature distribution deviates from 
that only slightly. 


Measuring Equipment 

The measuring instruments used were: Hypodermic-needle total-^head 

tube, Chromel— Alumel thermocouple, and hot-wire anemometer. 

The hot wires were nominally 0.00025— inch platinum etched from 
Wollaston wire. The etched platinum was soft solde ed to the tips of 
small steel needle supports. The anemometer heating circuits and 
amplifier were designed and built by Carl Thiele in 19^-1 • The circuit 
is so arranged that hot-wire time constants are determined by super- 
imposing equal alternating— current voltages at two frequencies upon the 
balanced direct— current bridge. 

The amplifier gain is constant to within ±2 percent over a frequency 
range from below 7 cycles to 7000 cycles. No check was possible at less 
than 7 cycles since that was the lower limit of the available oscillators 
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Compensation for hot-wire lag is achieved by a resistance— inductance 
network between two stages of the amplifier. For a normal range of wire 
time constants, the combination of hot wire and properly compensated 
amplifier is satisfactory over the entire flat range of the uncompensated 
amplifier. The amplifier output was read on an approximately critically 
damped Wall galvanometer, with a vacuum thermocouple. 

Mean total head (i.e., dynamic pressure in a free jet) and mean 
temperature were photographically recorded simultaneously, by means of 
an automatic traversing arrangement, used with the total-head tube and 
thermocouple . The total-head pressure line was run into a small copper 
bellows which tilted a mirror, thereby deflecting a narrow light beam 
upon a uniformly moving sheet of sensitized paper. The simultaneous 
recording of temperature on the paper utilized directly the light beam 
reflected from the mirrors of a sensitive galvanometer. The galvanometer 
was critically damped for all sensitivity settings. The possible errors 
arising from changes in "steady— state" conditions during a continuous 
unidirectional traverse were investigated by means of a few check runs 
in the opposite direction; these showed no appreciable difference. 

The automatic traversing was accomplished by means of a screw— driven 
carriage running horizontally along a steel track (fig. 2 ). The screw 
was rotated through a gear, worm, and belt drive by a reversible 
alternating-current motor with wide speed range, operating on a con- 
tinuously adjustable transformer. The sensitized— paper holder was 
mechanically connected to the moving carriage, so that the abcissas of 
the recording curves were equal to the true radial distance. Paper 
shrinkage in development was found to be negligible. Figure 4 is a 
typical record. The symmetry axes of total head and temperature are 
offset laterally because of the necessary distance between tube and 
thermocouple. 


PROCEDURES 

Mean Dynamic Pressure 


Mean dynamic pressure is effectively the total head in a free jet 
and was recorded photographically, as described in the previous section. 
The photo- record was faired and traced on cross-section paper. In 
order to reduce scatter, each traverse finally presented in this report 
has been averaged from three individual faired runs. Faired curves 
from typical repeated traverses are plotted together in figure 5 to 
illustrate the degree of disagreement. These particular runs were 
made with the old hot— jet unit. 

No correction has been made for the effect of velocity and density 
fluctuations upon the readings, but a brief analysis of this correction 
is given in appendix A. 
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Mean Temperature 

The temperature was also recorded photographically, as described in 
the previous section. Since the galvanometer deflection varied linearly 
with temperature, there was no reading error specifically attributable to 
the temperature fluctuations. 


Mean Velocity 

With the assumption of a perfect gas, velocity was computed from 
dynamic pressure and temperature. 


Axial Component of Turbulence 

In the cold jet, the axial component of turbulence was measured in 
the conventional manner with a single hot-wire, about 2 millimeters in 
length, set normal to the x— and r-directions. 

In the hot jet, ^r- was measured simultaneously with ^7- and -=r=> 
TJ e eu 

by the technique described in references 15 and 16. 


Throughout the report, no correction has been made for the error 
(in reading fluctuation levels) due to the fact that the fluctuations 
are not always small, as assumed in the theory of turbulence measurement 
with hot-wires. However, some extremely brief check measurements, 
reported in reference 12, indicate that the absolute values may be 
fairly good, even at levels above 50 percent. 


It is important to note that the mean velocity to which the hot 


wire responds is Up = |/u^ + V^ rather than simply U. Therefore the 
local turbulence level measured, and plotted in the figures, is 

actually =^-. However, this differs appreciably from ^ only at the 

% U 
outer edge of the jet, where the flow is not completely turbulent anyway. 

U f 

The distributions of — , however, are just that quantity at all points, 

Umax 


having been computed from 7^ — = =— 

Umax Umax 


u» U R 


The mean velocity distributions 


U 


plotted in the same figures as the fluctuation measurements are all — 


R 


^max 


and therefore, at the outer edge, are appreciably different from the 
curves obtained from total-head— - tube measurements and plotted in other 
figures. 
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Temperature Fluctuation Level 

Temperature fluctuation level was measured directly "by using a single 
wire as a pure resistance thermometer. 


Radial Component of Turbulence 

The radial component of turbulence was measured in the hot jet with 
an X— type meter made up of two wires (about 3 mm long), the voltage 
fluctuations of which were subtracted before amplification. Because of 
almost unavoidable unsyrametry, the meter was slightly sensitive to 
appropriate correction has been made. ^ 


•flu Correlation 


The flu correlation was measured simultaneously with 
as described in references 15 and 1 6 . 


U 


and 



uv and flv Correlations 

The uv and flv correlations were measured simultaneously, according 
to the general procedure given in the foregoing references. The extreme 
slowness of the measurements precluded the use of a simple inversion 
procedure for unsymmetry correction (i.e., repetition of readings at a 
point, with the meter rotated through l80°). Instead, the results were 
computed directly with the directional sensitivity calibrations of the 
individual wires. 


The mean squares of the fluctuation voltages of the two wires are 

°i 2 - * h 2 (ff + e i 2 (t) 2 + § + 2 °* 5 i f + 2 vi P 


e/ - <-4/ * r e *$f ♦ ^ ♦ a*r 8 § - a* f - f 


2 


NACA TN No. 1865 


11 


Since there are two unknowns, readings must "be taken at two different 
hot-wire sensitivities. Then the two pairs of differences give 



+ 2 (oqT! - -rjr + 2 ( a l s l + ^ 2 ) ^ + 2 ( 7 l 5 l 


+ ? 2 & 2 


\ uv 

V ~2 







+ 2 (a 1 , r 1 ' - VV) + 2 ( a l' 8 l’ + a 2* 5 2*) ^ 


+ 2 ( r l' S l' + ^2* 8 2*) ^ 


Thus, the u, v, and *6 sensitivities of each wire must be determined 
separately. From previous measurements, /—A , /^-\ , ^ and 

W \U/ 0U \U/ 

are known, and the foregoing pair of linear equations can he solved 
simultaneously for the two unknowns Ar~r and- 


eu 


U 2 


EXPERIMENTAL RESULTS 


The axial distributions of dynamic pressure, temperature, and 
velocity are given in figures 6 and 7 , for the constant-density jet 

Poo 

for an initial density ratio — ~ 2. Since all three variables 

p o 

decrease more rapidly when jet density is appreciably less than that 
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of the receiving medium, It is to he expected that the lowei^-density (hot) 
jet spreads at a wider angle. This is verified in figure 8, which shows 
the half total— head diameter 2r-^ and the half temperature diameter 2r* 

as functions of 2L, for these two cases. 


In figures 6 and 7 it can he seen that hoth dynamic and thermal 
potential cones are shorter for the low— density Jet. Since the annular 
mixing region outside of the potential cone is somewhat similar to a 
single two-dimensional mixing zone, the nature of the flow in this latter 
case, with flowing medium less dense than stationary medium, can he inferred. 

Figures 9 and. 10 show the dimensionless profiles of total head 
temperature, and velocity at a section 15 diameters from the orifice, 

P oo 

for — £ 1 and 2 . It is immediately evident that the shape of 
^o 

these functions is not appreciably changed at - = 15 by a doubling 

of the initial density ratio. The local maximum density ratios at 

P 

the section are ■=-£- « 1 and 1.3, respectively. A verification of 
Pmin 

the negligible change in total-head profile is given in figure 11, 
in which all the functions have been fitted at the point r = r-^. 

For a constant— density jet, it is well known that effective 
similarity of total— head (or velocity) profiles exists for a.1 1 values 

of £ greater than 7 or 8 (references 11 and 12 ). This is verified in 

figure 12 . Figure 13 shows that an initial density ratio of 2 is not 
sufficient to change this result appreciably. The deviation of the 

total— head distributions for ^ = 20 at large values of — is almost 

d 0 ^ 

certainly due to inaccuracy of the traced data at a section where even 
the maximum dynamic pressure gives only about a 1— inch deflection on the 
photo record. 

The fact that the density— ratio range covered in the present series 
of experiments is insufficient to lead to measurable deviation from 

f 

simple geometrical similarity is clearly shown by the constancy of — 

A . r l 

and — in figure 14 . However, again it should be pointed out that, 

although the initial density— ratio range is from 1 to 2, the range of 
maximum density ratio at this particular section is only from 1 to 1.3* 

Although r]_ has been used as a characteristic jet width in most 
of the previous figures, it is clear that a width definition of greater 
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physical significance can he made on the basis of momentum flow. The 
"momentum diameter" of the jet at any section 


A = 



is defined as the diameter of a Jet of rectangular density and velocity 
profiles (p = and U = n maY '), whose total rate of flow of axial 

momentum is the same as in the actual jet at that section. Since the 
momentum flow in a free jet is effectively constant, this characteristic 
diameter is particularly convenient. It is analogous to the momentum 
thickness of a boundary layer, except that the latter is based on 
momentum defect. 

The ratio of hot— jet momentum diameter at ^=15 to that of the 

constant-density jet at the same section (fig. 15 ) shows a definite 
increase in jet width as the jet fluid becomes less dense than the 
receiving medium. A few points obtained with the preliminary hot— jet 
equipment have been included because this simple unit was able to 
produce an appreciably higher initial temperature. The result of a 
rough theoretical analysis of the effect, presented in the next section, 
is drawn in the same figure. 

A comparison between the cold— jet axial u* distribution (fig. 1 6 ) 

and that for an initial temperature difference of 170° C (fig. 17) shows 
that the sharp rise in turbulence level comes a bit earlier for the low- 
density jet. Of course, this may be due to an inaccuracy in King*s 
equation and thus may not be a real effect; however, the difference is 
certainly in the direction that would be anticipated as a result of the 
more rapid spread and development of the mean motion pattern in the 
lower density case. 

The temperature fluctuation level as measured is everywhere lower 
than the velocity fluctuation level. 

Two qualitative aspects of the fluctuations were observed during 
these particular measurements: First, extremely 'regular velocity 

fluctuations were observed in the potential cone, just as illustrated 
in reference 13 ; second, the temperature fluctuations on the jet axis, 

roughly between j = 3 and 5 , were almost entirely one-sided, although 

the u fluctuations in that region apparently showed no such tendency. 
The reason for thisp behavior is not clear at the time of writing. 
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In the process of measuring 


u' 


in a flov with simultaneous 


U 


velocity and temperature fluctuations, the flu correlation results more 


flu 


Su y u > 


are 


or less as a byproduct; _ x i and the coefficient R 

e u 

o o 

plotted in figure l8. In spite of tremendous experimental scatter, which 
has even produced occasional negative values for #u, it appears 
that -9u > 0 in turbulent shear flow. The sign can be obtained by the 
reasoning (in rough analogy to the kinetic theory of gases) that uv 

^\t T 

has a sign opposite to that of ^ in a normal turbulent shear flow. 


Since $u is proportional to the turbulent heat transfer in the 
main flow direction, it is clearly not of as great interest as the iv 
correlation. It is, however, of use in the computation of the total- 
head— tube correction for the effect of velocity and temperature fluctua- 
tions (appendix A). No attempt has been made to draw a curve through 
the experimental points. 

A check measurement of the lateral u* distribution in the fully 
developed constant-density jet, in figure 19, shows reasonable agreement 
with earlier results (reference 12). 

Figure 20 contains lateral traverses of both u* and d * in the 

nr Q J -in ^ 

hot jet at -5- = 20. Again, — is everywhere slightly less than 

0 , U 

d* u 

but it should be noted that — > = in the outer part of the jet. 

e max ^max 

Perhaps associated with this relative behavior is the fact that, 

u* 

although — — has a definite local minimum on the jet axis, = is 

®max ^max 

flat in that region, within the accuracy of measurement. 

The coefficient R is apparently more or less constant over the 

flU 

central part of the jet (fig. 21 ). 


Considerably more complete measurement s have been made at — = 15 

in the hot jet. In addition to u f , d f , and du distributions 

(figs. 22, 23, and 24), there is also a traverse of the lateral component 

f yt 

of turbulence — and — . It may be noticed that, in agreement with 

U Umax 

figure 24 of reference 12, v 1 > u f in the immediate yicinity of the 
axis, and u f > v 1 elsewhere. 
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The uv and 3v correlations and correlation coefficients 


Ei 


■uv 


UV 

u*v* 


and E 

dv 


t3y 


are plotted in figure 25* There is a great deal of experimental scatter 
but the general behaviors are clear. Furthermore, it is obvious that 
K^ v > E uv , which is exactly what is to be anticipated from the fact that 

the rate of turbulent heat transfer is known (from the mean velocity and 
temperature curves) to be greater than the momentum transfer. 

Of course, because of the effect of density fluctuations, the 
turbulent shear and heat transfer cannot be expressed simply as — puv 
and -CpPflv, respectively. A brief calculation (appendix D, reference 16) 

leads to the approximate expressions 

__ _ U — 

t = -puv + p nr- vV 



In figure 2 6 , the shear stress and heat transfer are plotted in 
the dimensionless forms 
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ANALYTICAL CONSIDERATIONS 
Jet Width as a Function of .Initial Density Ratio 


As mentioned in the INTRODUCTION, integrated equations, along with 
dimensional reasoning in transforming the variables, can usually be made 
to yield solutions describing the integral behavior of a shear flow. 

The two-dimensional integrated momentum and mechanical— energy 
equations have been derived by Liepmann and Laufer (reference 6). An 
identical procedure leads to the axially symmetrical forms, omitting 
bars for mean motion. 

Momentum: 



pU^r dr — p b U b 2 b g + 


2 a U a 2a + M! 



(1) 


mechanical energy: 





1 

2 


Oa D a 3 a 


da 

— + 

dx 


1 

2 



b 

a 


= [ rU l]a ~l [ TT ^ ( 2 ) 

A similar treatment of the heat-transfer equation, with specific 
heat assumed independent of temperature, gives, again omitting the bars 
since all quantities are mean quantities, 



pURr dr - pJJJ b + pU„ 0 „a ^ + 
bbb dx K aa <3* 
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For the particular case of a free jet entering an infinite medium, 
at rest, these three equations simplify to 



dU 2 ! 


dr = Constant 


2 « 


w 


d_ 

dx 


pU 3 r 


dr = -2 



(5) 


pU0r dr = Constant = — 

2rt 


(6) 


where the integrals have been extended to infinity to permit the use of 
asymptotic velocity— and temperature— profile assumptions. If finite 
width functions are used, the integrals can he extended only to some 
judiciously chosen h. It should also he mentioned that in the integral 
continuity equation, which can he used to find the radial velocity Y^, 
the finite limit is essential. 

Equation (4) expresses the condition of constant total flux of 
momentum across all planes perpendicular to the jet axis in the absence 
of a pressure gradient. Equation ( 5 ) equates the rate of decrease of 
mean— flow kinetic energy per unit time to the rate of production of 
heat (in the laminar case) or of turbulent kinetic energy (in the 
turbulent case). 

Equation (6), expressing the constancy of heat flow across all planes 
perpendicular to the x— axis, is of course the thermal part of the energy 
equation, for the case in which heat production by viscous action is 
negligible. When the density variation is isothermal and due to the use 
of a jet fluid different from the receiving medium, equation (6) is 
replaced by what may be considered the conservation of mass— flow defect 
or excess: 



. . m' 

(p — Poo)Ur dr = Constant = — 


( 7 ) 
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and, in fact, equation ( 7 ) is easily obtainable from equation (6) with 
the equation of state of a perfect gas, assuming constant pressure . 

If there is introduced an appropriate dimensionless variable, along 
with reasonable assumptions for the dimensionless velocity and density 
profiles, equations (4), ( 5 ), and ( 7 ) can be solved for U^^ x), p m1ri (x), 

and the momentum diameter A(x). The integrated continuity equation can 
be used to solve for radial velocities. Although t can in principle be 
determined from the two assumed functions (see, for example, references 6 
and 12) it is sometimes more convenient to make a reasonable assumption 
for this distribution also. 


The purpose of the present analysis is only to obtain an approximate 

expression for — as a function of — , for a fixed value of it turns 
d P 0 fl- 

out that this can be done without the mechanical— energy equation. 

The velocity and density profiles are represented by dimensionless 
functions 


U 


II 


= f(q) and 


P ~ Pc 


m a.x 


P . - Poo 

min 00 


= s(t)) 


( 8 ) 


where q is the appropriate similarity parameter. It can be shown, by 
resulting inconsistencies in the two equations to be solved, that simple 
geometrical similarity does not exist. In fact, it seems physically- 
evident that the profile functions should vary with the local density 
ratio. Taking 



the form of F is obtained by substituting equation (8) and the definition 
of momentum diameter 

^ Pmin^inax = ^ 

Then, 

/ *\ p p A irA^ p 

vPmin — Poo/8j^max ^ ^2 4 A- 7 ! = ^ ^min^max 


into equation (4). 


Z oo 

C p “ + 
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which gives 


F = J2 


Ig + (!l - I 2 ) 


Pfeo 


-, 1/2 


p min 


where Ip and are pure numbers dependent upon the assumed velocity 

and density functions 


I-, = 


T]f 2 ( t \ ) dT) and 


I 2 = 


T|f 2 (T))g(T]) dTJ 


Thus, the proper similarity parameter is 




1/2 


I 2 + (!,_ - I 2 ) 


min 


(9) 


which is a function of x. If the problem were solved for p mlT1 (x) 

and A(x), t ) could then be expressed in its most suitable form, as a 
function of r, x, and the initial density ratio. 


It is of some interest to inspect the parameter t). For constant- 

density flow, r] = /2Ip and simple geometrical similarity exists. 

For heated turbulent jets it is found that Ip > I 2 , since g ( tj ) <_ 1.0 
everywhere. Thus, with p Tt11 - T1 appreciably less than , a particular 

value of t] is reached at smaller r (for fixed A) than in the 
constant-density case. 

Substitution of equation ( 9 ) into equation (4) merely gives the 
definition of A, which is now written with the explicit expression 
for M: 


p , U 2 A 2 = 
min max 


Po u o 


2 j2 


( 10 ) 


With equation ( 9 ) and the explicit expression for 


M, equation ( 7 ) becomes 


Pmin^inax^ C p min p oa) 

^ IpPco + Ip ( Pmin ~ Poo) 



( 11 ) 
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where 


X 3 


nf(n)g(Ti) dt) 


of 


Equations (10) and (ll) are to "be solved for 

P 0 Y 

at a particular value of Since there 

Poo d 


— as a function 
d 

are fewer equations 


than unknowns, some assumption must be made concerning the axial distribu- 
tion of either velocity or density. It is convenient to make the 
simplifying assumption that the dimensionless axial density distribu- 


tion 


/ p min p oo\ 

V p o - P «J 


is independent of the initial density ratio. 


U 


Then 


max 


from equation (10) is substituted into equation (ll). 


and all p^.^s are expressed in terms of 


' Pmin ~ P K 

V P 0 “ P oo, 


After further 


approximation by the use of the first two terms of binomial, expansions, 
the following result is obtained: 



where is the momentum width at the same axial distance for the 

constant— density jet, and where the desired — is obtained by choosing 
' p min “ P~ x d 


according to temperature measurements in the constant-density 


case. Thus it is clear that, as may he intuitively anticipated, the jet 
spreads at a wider angle when the jet fluid is less dense than the 
receiving medium, and vice versa. 
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In 

A numerical value for ^ is obtained by computation from velocity 
and temperature distributions in the measurements of the constant— density 


case. These give = 0.727 and equation (12a) is plotted in figure 15 

*1 

for comparison with the experimental results. Because of the rather 
severe approximation, agreement is not good. The complete solution of 
equations (4 ) , ( 5 ), and (j), however, is rather laborious. 


Measure of Turbulent Prandtl Number 

One of the more important and interesting problems in the study of 
turbulent shear flow is the relation between .momentum transfer and heat 
transfer. 

For laminar flows, the ratio of momentum to heat transfer is directly 
computable from the kinetic theory of gases and is conventionally 
described by the well-known Prandtl number 


a 


k 


For nearly perfect gases, o is only slowly variable with temperature 
over a wide range since Cp is nearly constant and u and k vary 

nearly proportionally. For effectively isothermal laminar flow, Cp, n, 
and k are individually constant over the entire flow field. 

The quantity analogous to a, for turbulent flow, is 



( 13 ) 


which may be termed the "turbulent Prandtl number." The "eddy viscosity" 
and "eddy conductivity," as defined by Bryden (reference 19 ) are 


— dU 
-puv = 

dr 


— „ de 


( 14 ) 


lln this section all U' s and 9 ' s are mean values. 
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and if the density fluctuations are negligible, as is the case within 
the accuracy of the present experiments. 



Q = 



Since e and (B are themselves rather artificial quantities, introduced 
merely in overstrained analogy to laminar flow, it is not to he expected 
that C-j, has any fundamental physical basis, other than the exact 
definition given in equation (15 ) . However, in the absence of a satis- 
factory analysis of turbulent shear flow, it serves as an empirical 
measure of the ratio of momentum transfer to heat transfer. Although e 
and p vary from point to point, even in an isothermal turbulent shear 
flow, the variations in their ratio may be of distinct interest. 

An excellent discussion from the point of view of the mixing— length 
theories has been given by Dryden (reference 20 ) and need not be repeated 
herein. 

From equations ( 13 ) and (l^), the local turbulent Prandtl number in 
terms of measurable quantities may be written as 


(15) 


Se 

uv Sr 
3v SU 
Sr 


or, in more convenient computational form. 



(15a) 


On the basis of the experimental results presented in the previous 
section, two different measures of can he computed: First, the 

effective turbulent Prandtl number for an entire section of the jet 

d, can be computed from a comparison between the width ratio of mean 
"G 


NACA TN No. 1865 


23 


velocity and. temperature profiles in the turbulent jet and the same 
ratio as computed theoretically for a laminar jet. Second, the direct 


measurement of 



and 


dv 


\,9 U 

\ max maxy 


permits direct computation 


course, this 


of the 0 . distribution, according to equation (15a). Of 


latter result will not be very accurate, both because of the uncertainty 
of King*s equation and because of the excessive experimental scatter. 


Effective turbulent Prandtl number .— Effective turbulent Prandtl 
number provides another illustration of the usefulness of the integral 
equations when a complete solution is either not necessary or not feasible. 
For simplicity, the analysis is restricted to the constant-density case. 


In addition to equations (4), (5), and (6), for the round jet an 
equation analogous in form to the mechanical— energy equation (2) can be 
derived by multiplying the heat— transfer equation by U before integration 



r - / dr = - 



Sr p Jo 


rr 


Se 

Sr 


dr 


(16) 


The physical significance of equation (l6) is not so clear as that of 
equation ( 5 ). However, it may be regarded as a device for retaining Q 
(which disappears in the integration of the heat^ transfer equation) . 

The solution for the velocity field of a laminar round jet has been 
given by Schlichting (reference 20 ) and by Bickley (reference 21 ) , but 
the temperature field apparently has not been treated. For laminar flow. 


SU 

t = u — and 

Sr 


are substituted into equations (5) and 


Q = k 
( 16 ), 


Se 

Sr 

giving 


d_ 

dx 



TPr dr = — 2v 



dr 


(IT) 
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and 


+ fr Tr * (18) 

It is to be expected that a good measure of the Prandtl number will 
be the difference in momentum and thermal diameters of the jet. The latter 
is defined so that 


pU0r dr = N 


Then equations (4) and (6) are replaced by the definition of A and 
the foregoing equation with p = Constant, that is. 


nAe 


(3 min^max ®max ^ 


it 2.2 
^max ^ 


4M 

itP 


(19) 


^max^maxA©^ - 


4n 

jtpc T 


( 20 ) 


With constant density, simple geometrical, similarity exists, and there 


may be introduced the dimensionless variable q* = along with the 


functional assumptions 


U 


U. 


= f ( q * ) and 


max 


6. 


= h(q* ) . 


max 


Substitution of q* and f into equation (17), followed by 
substitution from equation (19) into the result, gives the differential 
equation 


dU, 


max 


= -2jt 


E 2 ii 


XL 


I, M max 


dx 


( 21 ) 
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Substitution of q*, f, and g into equation (l8), followed "by 
substitution from equation (20 ) into the result, gives the differential 
equation 


The I f s are dimensionless integrals of f, g, and q*. Of course, 
equations (2l) and (22) can he solved for UmaxW ani e mpY (x), whence A(x) 

and Aq ( x) are obtainable from equations (19) and (20). However, these 
results are not of primary interest herein. 

Equating equation (21 ) to equation (22) gives 


Without evaluating the I*s individually, the product is obtained 
by considering as boundary condition the solution for 0=1, which 
means that momentum and heat— transfer rates are identical. In this 
limiting case Ag = A, since the velocity and temperature distributions 



( 22 ) 



(23) 


Division of equation (20) by equation (19) gives 



which, with equation (23), yields 



l 1 / 2 


(2b) 


must be identical. Therefore, 




( 25 ) 
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or, if this method la used to dete rmin e 


a instead of to predict 


Ae 

A' 



(25a) 


An identical treatment of a constant— density turbulent jet on the 
crude assumption of constant e and p results , of course^ in 



(26) 


where 0^ is the over-all effective Prandtl number for the turbulent— Jet 
cross section. 

For air, in the normal range of room temperature, the considerable 
scatter in thermal conductivities obtained in the various experimental 
determinations leads to an approximate range of possible true (laminar) 
Prandtl numbers, 0.70 < 0 <0.76. A tabular and graphical summary of 
"best" measured values up to 1938 has been made by Tribus_and Boelter 
(reference 23):_ They give, for example, 0 = 0.708 for T a = 302° abs. 
and 0.701 for T a = 3^8° abs. These particular values have been chosen 

for comparison with two computations of the turbulent measurements . 

For the present turbulent jet, at ^-=15 and initial temperature 

differences of 0 = 15° C and 300° C, the width ratios are — = 1.102 

_ A 

and 1.107, respectively. From equation (26), these give 0^. = 0.701 
and O.69O. With the two orifice temperatures given, the average absolute 
temperatures across the section are about 302° abs. and 3^8° abs. (these 

numbers are actually T r + jy 0 mBY ). This gives the extremely interesting 

result that the average of effective turbulent Prandtl number in a round 
jet is equal to the laminar value, within the probable accuracy of the 
present measurements. 


Turbulent Prandtl number distribution.— Faired curves for the 


uv 


U ^ 
u max 


and — 


■d v 


0 U 
max max 


measurements of 


(fig. 25) have been used, along 
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with the dimensionless velocity and temperature profiles, to compute 
local values of cr-^ in accordance with equation (l5a) # The values have 
"been plotted in figure 27 ., along with a heavy line giving 0 ^ as 
A n 

computed from Near the axis and near the jet edge, the computation 

becomes more or less indeterminate. 

The agreement between a^_ and the average value of the points 

is extremely good. In fact, it is rather better than was anticipated, 
because of the considerable scatter of the fluctuation measurements and 
the possible inexactness of King's equation. 


DISCUSSION 

Errors 


Most of the sources of error and the corrections, both applied and 
not applied, have been covered in the previous sections. However, it 
may be well to outline and extend the previous remarks. 

Total-head— tube and thermocouple data .— The total-head— tube readings 
have not been corrected for the effect of velocity and density fluctua- 
tions (appendix A) . Also, since the range of possible sensitivities was 
quite limited and since the response is proportional to U^, the deflec- 
tions near the jet edge were small and difficult to obtain accurately. 

A certain amount of experimental_scatter in al ~l measured quantities 
resulted from small fluctuations in U Q , due primarily to the partly 
stalled condition of the blower, which was designed to operate against a 
somewhat smaller pressure rise. 


— T i 

The temperature distributions for small 8 Q and at large and — 

were naturally susceptible to considerable scatter due to slow temperature 
fluctuations of perhaps a fraction of a degree in the laboratory. 

Hot-^wlre measurements .— The usual sources of error (see discussions 
by Dryden (reference 23) and. Simmons (reference 24)) for hot-wire 
measurements were encountered in the present investigation. 

Two additional difficulties were attributable to the particular flow 
configuration: First, the extremely high turbulence level present in 

a free jet is certainly outside of the small disturbance assumptions 
basic to the standard hot-wire response calculations. However, some 
rather brief measurements (reference 12) indicate that opposite errors 
may reduce this error considerably. Second, the extremely large degree 
of fluctuation introduces some uncertainty in the bridge balancing for 
the setting of hot-wire resistance. 
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The principal additional source of error inherent in the hot— jet 
measurements was the extreme difficulty of determining R a , the unheated 
wire resistance at local ambient temperature T a . Since the d* 

sensitivity of a wire is proportional to R a — R r (where R r is the 

unheated resistance at T r ), it is clear that, in regions of relatively 


small 0, a small error in R c 


it turns out, also in 

0U 


may cause a large error in 4L-, and as 

0 


The complexity of the computational procedures 


in these new applications of the hot-wire anemometer, often involving 
small differences between two large quantities, each with a fairly large 
amount of experimental scatter, has obviously magnified this scatter 
immensely in' the final results. 


Because of the foregoing factors, it was felt that hot-wire length 
corrections were not appropriate; the scale of the turbulence was of the 
order of 0.5 to 1 centimeter, and the scale of temperature fluctuations 
was probably of the same order of magnitude. 


However, it should be emphasized that free jet measurements represent 
the most stringent kind of test for the new measuring technique. It is 
certain that similar investigations in channels and boundary layers, for 
example, will show more consistent experimental points. 


Variation of Jet Spread with Density 


As pointed out, the considerable divergence between computed and 


measured curves of A. against ^ (fig. 15) is probably due to the 

A 1 p o 

very rough nature of the simplifying assumption used. Apparently the 
curve shapes agree reasonably well. This theoretical curve is computed 
particularly for the case in which density difference is obtained by 
heating; the relation between heat and momentum transfer determines the 
relation between the velocity function f(q) and the density function g(T)) 


and thus determines the ratio 


*i 


in equation (12). In cases where the 


density difference is obtained by the use of different gases, may be 

A p X 1 

different, so that the curve of — against — will not be the same. 

A 1 p o 

However, it does not seem likely that material diffusion in turbulent 
shear flow is much different from momentum and heat diffusion, 2 so the 
final result should be fairly close to that computed herein. 


p 

At the Seventh International Congress for Applied Mechanics, 

B. G-. van der Hegge Zijnen and J. 0. Hinze reported careful measurements 
showing heat and material diffusion to be effectively the same in a 
round turbulent jet . 
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Fluctuation Levels 


There are two characteristics of the measured fluctuation levels 

u * t3 ® 

that are worthy of note. The first is that ■= > — everywhere in the 

U Q 

jet. Although there must he some uncertainty since the complete physical 
form of King f s equation has never been verified (appendix A, reference 16 ),, 

u f 

the fact that — in the heated jet is measured at about the same value 
U 

as in the cold jet seems evidence that the results are reasonable. An 

-6 u* 

additional bit of evidence is the result that — > — — 


part of the jet 


the axis where 




max 


II 


max 


whereas 


in the outer 
u* _ 


> 


near 


U, 


max 


ft 


max 


This is quite consistent with 


the result arrived at "by the conventional qualitative phenomenological 
approach which considers a u fluctuation (or a i3 fluctuation) arising 
from a v fluctuation with a fixed mixing length. Thus, if this (purely 
fictitious, hut sometimes convenient) mixing length were about the same 
for velocity and temperature fluctuations, u* would be relatively greater 
than -O’ where the U profile was relatively steeper than the 0 profile, 
and vice-versa. Of course, this does not necessarily account for the 


slight local minimum in 


d' 


on the jet axis. 


max 


Relative Transfer of Heat and Momentum 


The a ^ distribution shown in figure 27, although agreeing in 


average with cr, computed directly from 


*0 


does not necessarily prove 
Both the 


t 'A 

any significant result with regard to the function a^-(r), 
experimental scatter in the uv and t9v measurements and the difficulty 
of making accurate slope determinations from the mean velocity and 
temperature profiles permit only the conclusion that a ^ is roughly 
constant across the main part of a round turbulent jet. 


That the effective turbulent Prandtl number might possibly be the same 
as the laminar Prandtl number was apparently first suggested by Dryden 
(reference 19), though not with any particular physical justifica- 
tion. In fact, it appears that no one has yet put forth any rational 
hypothesis , e.ither rigorous or phenomenological, to predict the relative 
rates of heat and momentum transfer in turbulent shear flow. 
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Since the vorticity transport theory is apparently the only analytical 
approach that has led even to the qualitatively correct result, it seems 
probable that simple momentum-exchange considerations will be inadequate. 


SUMMARY OF RESULTS 


From measurements of mean total-head and temperature fields in a 
round turbulent jet with various initial temperatures, the following 
enumerated results can be stated with confidence: 

1 . The rate of spread of a round turbulent jet increases with a 
decrease in the density of jet fluid relative to receiving medium. When 
the density difference is obtained by heating, the effect is as shown by 
the experimental points in figure 15. 

2 . Up to a local maximum density ratio Poo/p ffllT1 of about 1.3 

simple geometrical similarity still exists in a fully developed jet, 
within the accuracy of the present measurements. The total-head and 
temperature profile functions are effectively the same as in the constant- 
density jet. 

3 . The average or effective "turbulent Prandtl number," for a 
section in the fully developed jet, is very nearly equal to the laminar 
Prandtl number. 

4 . The local turbulent Prandtl number at a Section in the fully 
developed jet, away from axis and edge, is roughly constant. In addition, 

d* u* 

with somewhat less certainty, it has been found that — < — every— 

9 t U » 

where in the jet and that the relative magnitude of ~ — and ^ — 

0ma * V 

varies roughly as the relative magnitude of the slopes of dimensionless 
temperature and velocity profiles. 


California Institute of Technology 

Pasadena, Calif., August 18, 1947 
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APPENDIX A 

EFFECT OF SIMULTANEOUS VELOCITY AND DENSITY FLUCTUATIONS 
ON READING OF A TOTAL-READ TUBE 


Neglecting Doth static— pressure fluctuations and the effect of 
lateral velocity fluctuations, the tube reading is 


P = | (p + p* )(U + u) 2 


(Al) 


Therefore, since pUu = 0, 


P = i- p (u 2 + ^ + Up*u + p*u 2 (A2) 

For small fluctuations the last term can be neglected. Also, since the 
density fluctuations of interest herein are due to temperature fluctuations. 



is substituted into equation (A2), which finally gives for the true mean 
dynamic pressure in terms of the tube reading: 



P 


1 + 


U 2 



(A3) 


It is interesting to note that, since -flu is apparently ordinarily 
positive, the errors due to the velocity and the temperature fluctuations 
.are in opposite directions. 

The justification for neglecting the effect of lateral fluctuations 
in this approximate calculation is that its effect on the magnitude of 
the instantaneous velocity vector is at least partially balanced by the 
directional sensitivity of the tube . 
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(b) Close-up of orifice and traverse track. 
Figure 2 .- Test setup. 


(a) Complete hot -jet unit. 
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e 0 = 15° C 




r, in. r,in. r,in. r, in. 

Figure 3.- Initial temperature and total -head distributions. 1-inch round jet. q 0 = 40 millimeters of 

alcohol. 
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Figure 4.- Sample photographic record of mean total head and temperature. 1 -inch heated jet. Lateral 

offset of records is due to instrument spacing. 
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Figure 5.- Experimental scatter. 1 -inch heated jet. x = 15 inches; e Q = 280° C. 
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Figure 6.- Axial total -head distributions. 1-inch round jet. 
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Figure 7.- Axial velocity and temperature distributions. 1-inch round jet. 
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Figure 8.- Spread of heat and momentum. 1-inch round jet. 
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Figure 10.- Temperature and velocity distributions. 1-inch round jet. 
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Figure 11.- Check on similarity of total -head distribution. 1-inch hot jet. j = 15; 
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Figure 12.- Check on similarity of total -head distributions. Constant -density jet. e Q = 15° C. r-^ = r 

at which q = |q max . 
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Figure 13.- Check on similarity of total -head distributions. 1-inch hot jet. 
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Figure 14.- Comparison of characteristic lengths. 1 -inch heated jet. ^=15. 
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Figure 15.- Variation of momentum diameter with initial temperature or 

density ratio. Hot jet. ^ = 15. 
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Figure 16.- Axial velocity fluctuation levels. 1 -inch round jet. 9 0 = 15° C. 
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Figure 17.- Axial distribution of fluctuations. 1-inch hot jet. 
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Figure 18.- Axial distribution of tfu correlation. Hot jet. 0 Q = 170° C. 
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Figure 19.- Check measurement of lateral u' distribution. Cold jet. • | = 20. 


No. 1865 


VJ1 



NACA TN No . 1865 


NACA TN No . 1865 


55 


r / r 2 



Figure 21 . ~ The correlation. Hot jet. 2.= 20; e 0 = 170° C. 
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Figure 23.- Hot -jet fluctuation levels 
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Figure 24.- Hot-jet 4u correlations. 4 = 15; e 0 * 170° C. 
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Figure 25.- Hot -jet double correlations, j = 15; 0 O 


170° C. 


r / r 2 S' 



| - 15; e 0 - 170° C. 


Figure 26.- Hot -jet momentum and heat transfer. 
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Figure 27.- “Turbulent Prandtl number.” 1 -inch hot jet. j = 15; 

0Q = 170° C. 
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